Complete temporal characterisation of asymmetric pulse compression in a laser wakefield 
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We present complete experimental characterisation of the temporal shape of an intense ultrashort 200-TW 
laser pulse driving a laser wakefield. The phase of the pulse was uniquely measured using (second order) 
frequency resolved optical gating (FROG). The pulses are asymmetrically compressed, and exhibit a positive 
chirp consistent with the expected asymmetric self-phase modulation due to photon acceleration/deceleration in 
a relativistic plasma wave. The measured pulse duration decreases linearly with increasing length and density of 
the plasma, in quantitative agreement with the intensity dependent group velocity variation in the plasma wave. 



PACS numbers: 52.38.Kd, 41.75Jv, 52.50 Jm, 52.65.Rr 

Laser wakefield accelerators have now demonstrated the 
production of quasimonoenergetic electron bunches to GeV- 
scale energies over only ~ cm lengths (TJEl. Though remark- 
ably short for an accelerator, these distances are still many 
times the Rayleigh length, zr, of the driving high-intensity 
laser pulse, so that the non-linear interplay between laser and 
plasma becomes important. This leads to a wealth of in- 
teresting phenomena such as self-focusing OIH, self-phase 
modulation |H or photon acceleration [ 7 ] and pulse short- 
ening O |3. In combination, these processes can result in 
a 2 (intensity) amplification of the laser pulse El [8). Here 
a = eE/(m e uoc) is the normalised vector potential for a 
laser with electric field E and angular frequency c^o- Laser 
wakefields are driven by the ponderomotive force of the laser, 
F p = — ^m e c 2 Va 2 , where 7 2 = 1 + a 2 / 2 is the relativis- 
tic factor due to the transverse quiver of electrons in the laser 
field. Hence, determining the evolution of a 2 is of vital impor- 
tance in understanding wakefield accelerators. For example, 
a 2 amplification was essential in the first demonstrations of 
monoenergetic beam production ifTOl . 

For the interaction to extend to many Rayleigh lengths, the 
laser must be prevented from diverging naturally. This can 
be through the action of an external guiding structure [lj, or 
simply through a combination of ponderomotive and relativis- 
tic self-focusing Q. For sufficiently intense ultrashort laser 
pulses (a > 1), the laser will quickly focus to a matched 
spot radius, w rn « ^^/ao A p where X p = 2ttc/uj p , uj p is 
the plasma frequency and ao is the peak vacuum value of a 
l4lfm. Once matched, transverse variations in spot size play 
little role in increasing a 2 . 

However, as the laser sits in the density variation created 
by the plasma wave, it experiences a spatiotemporally vary- 
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ing refractive index rj = (l — cj 2 /(7cjq)) . For the front 
of the plasma wave, the density increase as plasma electrons 
are pushed forward is mostly compensated by the intensity de- 
pendence of t\ (through the dependence on 7) [12]. However, 
beyond this, both density decrease and further increase of 7 
cause r] to increase. This results in an increase in group ve- 
locity v 9 « 77c, and thus compression of the pulse increasing 
towards the back of the plasma wave. 



Associated with pulse compression must be an increase in 
spectral bandwidth. At the front of the pulse, density and 
intensity variations cancel to limit the amount of frequency 
shifts. However beyond the first maximum of the plasma 
wave, both decreasing density and increasing intensity cause 
an increasing 77 and thus a reciprocal decrease in phase veloc- 
ity v<p. Hence, the majority of the pulse can be redshifted. For 
longer pulses (cr « X p ), both increasing plasma density and 
relativistic self-phase modulation (decreasing 7) act together 
to cause strong blueshifts at the back of the pulse. 

At high intensity and including multidimensional effects, 
this picture is complicated further. The plasma wave becomes 
non-linearly steepened, and the centre of the plasma wave 
evolves into a caviton-like structure with relatively small vari- 
ations in 77 within it Hl3l . As a result, laser energy mostly 
moves forward towards the front of the caviton, where it is 
rapidly redshifted and thus slips quickly back within the wave 
frame. This can cause pulse steepening at the front of the 
pulse, and reduced steepening at the rear. 

In any case, the outcome of these complex dynamics is gen- 
erally compression of the laser pulse and possible increase in 
a 2 . Along with the exact pulse shape, this will determine F p 
and thus wakefield growth. Pulse compression of laser pulses 
in a laser wakefield has been reported previously [ 9 ]. Com- 
pression of a 30 fs pulse to ~10 fs, was measured with a 
second-order autocorrelator. However, this method does not 
give phase information and thus does not allow a complete 
diagnosis of the pulse shape. 

In this letter we report on the complete temporal (amplitude 
and phase) characterisation of short (r ~ 45 fs), relativisti- 
cally intense (/ « 2 x 10 19 Wcm -2 ) laser pulses after interac- 
tion with dilute plasmas such that the initial pulse was shorter 
than the excited plasma wave wavelength (cr < X p ). Spectral 
broadening, photon acceleration and asymmetric pulse short- 
ening have been observed. The dependence of these nonlinear 
effects on plasma density and interaction length have been in- 
vestigated, detailing the first parametric investigation of tem- 
poral pulse evolution in a laser wakefield. 

The experiment was conducted with the Astra- Gemini laser 
at the Rutherford Appleton Laboratory delivering laser pulses 
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at central wavelength Xl = 800 nm with energies up to 12 J 
in a FWHM duration of r = 45 fs. The 200 TW laser pulses 
were focused by a f/20 off-axis parabolic mirror to a spot size 
of (wo = 22.0 ± 0.6) fim FWHM, with corresponding confo- 
cal parameter of zr — 1 mm. For these parameters, ao « 3.0 
in vacuum. The interaction with a helium gas jet with elec- 
tron densities n e = 1.2 — 6.6 x 10 18 cm -3 and three different 
lengths I = 4, 6 and 8.5 mm were investigated. n e and I were 
determined by interferometry with a transverse optical probe 
El . Under these conditions the intensity modification was not 
sufficient to produce electron beams via self-injection. 

The transmitted laser pulses were collimated by a large 
aperture (//10) spherical mirror and imaged onto a 12-bit 
CCD camera to confirm the quality of self-guiding. The full 
beam was also focused onto the slit of an optical spectrome- 
ter sensitive over 300-1000 nm, that was absolutely calibrated 
using a white light source. An absolutely calibrated photodi- 
ode measured the transmitted laser energy. The temporal de- 
pendence of intensity and phase were measured using second- 
order frequency resolved optical gating (FROG) based on sec- 
ond harmonic generation (SHG) by guiding a small part of the 
transmitted laser (close to the centre of the beam) into a SHG- 
FROG device (GRENOUILLE by Swamp optics 1 14]) with 
12 fs temporal and 4nm spectral resolution. The amount of 
glass in the beam, such as vacuum windows, reflective neutral 
density filters and the long focal length lens for recollimation, 
was kept to a minimum and did not exceed 4.5 mm in total. 

In SHG-FROG, a second-order autocorrelation is spectrally 
dispersed in the direction vertical to the time axis. The result- 
ing 2D map - the FROG trace - thus contains both amplitude 
and phase information of the laser pulse. This allows the pulse 
shape to be uniquely determined, which is impossible from 
the autocorrelation alone. The pulse shape is established us- 
ing a multiparametric fitting algorithm, referred to as retrieval 
(FROG3 lfT6l ). In the retrieval, the temporal axis is discretised 
and, starting from an initial guess, the amplitude and phase of 
the laser at each point are varied until the measured FROG 
trace is best fitted by the synthetic one. The error is defined as 
the square root of the sum of the squared differences in each 
pixel of the measured and synthetic FROG. 

Experimentally obtained SHG-FROG traces can be asym- 
metric with respect to the time axis due to unwanted spatial 
chirp and pulse front tilt. Although we did not observe appre- 
ciable asymmetries, the FROG traces were symmetrised prior 
to the retrieval in order to improve the contrast and access a 
larger data set. The analysis was restricted to retrievals with 
an error smaller than 2% and for which the total spectrum 
measured in the spectrometer was in good agreement with the 
retrieved spectrum. This gives confidence that the FROG re- 
sults are representative of the complete pulse and that no lim- 
itation due to spectral clipping occurred. Finally, fields were 
corrected for propagation through glass, vacuum window and 
filters in the beam path. The pulse was attenuated by reflec- 
tion off glass wedges before passing through any material so 
that only linear dispersion in the optics need be accounted for. 

The complex laser fields are most usefully displayed in 



the form of a Wigner-distribution W(t, uj) fT31 which rep- 
resents the distribution of photons in (u,t) phase space. 
The marginals of the Wigner distribution, j W(t, uj) du and 
fW(t, uj) dt, give the temporal intensity and spectrum of the 
pulse and its first normalised moment gives the instantaneous 
frequency. Therefore, in a single 2D map it is possible to vi- 
sualise the most relevant information about the pulse. 

The remaining uncertainty of the direction of time due to 
the symmetric nature of the SHG-FROG was removed using 
physical considerations. First, a FROG trace of the pulse was 
obtained without gas, but with an additional 4 mm of glass 
in its path (Fig.[T^). The artificially introduced positive chirp 
(increasing frequency with time) increases the total chirp mea- 
sured (including the effect of the fixed glass in the system). 
Fig.[lJ shows the corresponding Wigner distribution after cor- 
recting for the glass in the beam path. We obtain an essentially 
bandwidth limited reference pulse with its temporal profile 
shown in black. The wings in the temporal profile are a com- 
mon artefact in high power ultrashort laser systems and has 
little influence on the physical processes to be described. 

Any similar positive dispersion due to the interaction would 
also increase the total chirp. Negative dispersion however 
would counter the positive chirp inherent in the system. So 
despite leading to the same pulse duration after the interac- 
tion, this would result in a quite different total chirp in the 
diagnostic. By varying the experiment parameters in small 
steps, so that the transition from the previous conditions were 
smooth, the change in the chirp and therefore the direction of 
time could be determined for each retrieved pulse. 

The pulse was measured after propagation through a plasma 
with relatively constant density n e G 1.2. .2.3 x 10 18 cm -3 , 
low enough that the laser (cr ~ 13.5 fim) did not extend be- 
yond the first plasma wave period (X p G 30.. 22 /im). For the 
shortest propagation distance I — 4 mm (Fig. []J>,f) we ob- 
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FIG. 1: (color online) Length scan at low electron densities n e G 
(1.2. .2. 3) x 10 18 cm" 3 . FROG-traces (first row) and Wigner- 
distributions of retrieved and corrected fields (second row) for / = 
0,4,6 and 8.5 mm from left to right. Overlaid black lines correspond 
to temporal intensity distribution, t < represents front of the pulse. 
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FIG. 2: (color online) Pulse profiles for n e — (a) 2.3, (b) 3.8, (c) 
4.7, and (d) 6.6 x 10 18 cm~ 3 for I = 4 mm. The black vertical lines 
mark the plasma period (2ir/cj p ), i.e. expected positions of electron 
density maxima for the respective densities. 



serve a positive chirp as expected. Noticeable though is the 
faster steepening at the back of the pulse (at positive times). 
The front (negative times) stays comparable to the initial pulse 
and the total pulse-length is reduced. For a longer distance of 
I = 6 mm (c,g), the front is depleted and for the longest dis- 
tance I = 8 mm (d,h) the pulse becomes near-transform lim- 
ited with a markedly reduced pulse duration (r = 27 ± 3 fs). 

Similar compression is observed for increasing n e for fixed 
propagation distance. Fig. [2] shows the Wigner-distributions 
and pulse profiles for the shortest interaction length (I = 
4 mm), but this time for electron density increasing from 2.3 
to 6.6 x 10 18 cm- 3 (X p e 22..13/im). The lowest density 
n e = 2.3 x 10 18 cm -3 (Fig. [2^) corresponds to the situation 
described in Fig. [lj>,d. With increasing n e , the pulses show 
increased bandwidth and shorter pulse duration. In contrast to 
the length scan at low density, the red- shift at the leading edge 
of the pulse is now more prominent. At n e = 6.6 x 10 18 cm -3 , 
the pulse is again almost transform limited with a broadened 
spectrum that supports the pulse duration r = 18 ± 3 fs. The 
pulses also exhibit enhanced structure with increasing density, 
with additional pre- and post-pulses observed. As a guide, the 
duration of one plasma period is indicated by the vertical lines 
at the top of each panel in fig. [2] The appearance of multiple 
pulses and pulselets indicates the modulation of wings of the 
initial pulse by the plasma wave. 

Fig. [3^, b quantify the pulse shortening. For constant den- 
sity n e ~ 1.8 x 10 18 cm -3 the FWHM duration reduces lin- 
early from r = 44 fs to 27 fs when the plasma length is step- 
wise increased to I = 8.5 mm. The rate of compression is 
« 2.0 ± 0.7fs/mm. Similarly, pulses are shortened from 
44 fs to 18 fs for a constant I = 4 mm and n e increased from 
2.3 to 6.6 x 10 18 cm -3 giving a constant compression rate 
^4.4±i; 8 fs/(10 18 cm- 3 ). 

Fig.[3j:,d show the variation in energy transmission, El ob- 
tained from two complementary measurements; the integrated 
spectrum and the absolutely calibrated diode. El decreases, 
with increasing propagation length and electron density. De- 
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FIG. 3: (color online) Dependence of pulse duration r on (a) prop- 
agation length I for n e = (1.2.. 2. 3) x 10 18 cm -3 and (b) electron 
density n e for I = 4 mm. Best fit lines give compression rates of 
2.0fs/mm and 4.4fs/(10 18 cm -3 ) respectively, (c), (d) shows trans- 
mission in terms of laser energy (full circles) and peak power (open). 



noting the normalised temporal pulse profile as f(l,n e ,i) 
(I = n e = being the reference pulse), the transmitted en- 
ergy, £x(Z, n e ) = J P p (l, n e )f (Z, n e ? t) dt, can then be used 
to calculate the relative peak power 
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where T(Z, n e ) = El(1, n e )/Ex(0, 0) is the energy transmis- 
sion. Fig.|5j:,d show a substantial initial drop in P p to « 30 % 
of its original value for I = 4 mm. This can be attributed 
to the initial self-focusing, which traps only a fraction of the 
original laser energy in the matched guided pulse, especially 
for non-ideal beams |4 |. However for further propagation, the 
peak power changes only slowly over the range of parameters 
investigated. Evidently, under these conditions, pulse com- 
pression compensates for the energy loss of the pulse in con- 
tinuously driving the wakefield. 

The pulse compression can be understood by considering 
two points in the pulse separated by cr having different group 
velocities v g \ and v g 2- Their separation therefore changes af- 
ter propagating in z direction over a distance Al according 
to cAr/Al (v g2 - v gl )/c (r/ 2 - Vi) ~ cr o§i> wnere 
ro is the initial pulse duration. As we will show, as a rea- 
sonable first approximation, ^ can be considered to be con- 
stant over the part of plasma period in which the majority of 
the laser energy sits for a pulse of duration cr w X p . Hence 
|f ~ (Vmin ~ Vmax)/cr. For underdense plasma n <C n c , 
we have 77 « 1 — \{n e /^n c ), where n c = meQUJ 2 /e 2 is the 
critical density. Hence for sufficiently intense pulses; a > 1; 
Vmax ~ 1; Vmin - 1 - (n e0 /2n c ). This gives a variation of 
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pulse duration on initial density n e o and I, r = tq - 
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the two scans performed, we find compression rates At / Al = 
—n e o/2cn c = — 1.7fs/mm for n e o = 1.8 x 10 18 cm -3 , and 



Ar/An = -l/2cn c = -3.8fs/10 18 cm- 3 with I = 4mm, 
in good agreement with the measured rates shown in Fig. [3] 

The compression of a laser pulse in a laser wakefield has 
been studied through simulation in 1 and 2D for similar pa- 
rameters by Gordon et al 0. These simulations display 
many of the same traits as our observations, in particular, re- 
producing a positive chirp for short propagation distance at 
low density. However to gain better insight into the pulse 
evolution, we consider the 1 -dimensional quasistatic wave- 
equation for a pulse driving a wakefield |12|. Writing the 
normalised vector potential ol in terms of a slowly varying 
envelope and a carrier frequency term fc = ^o/c, t) = 
^a(^r) exp(zfcoO + c.c, the wave equation for the complex 
laser envelope reads 
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where ^ — z — ct and r = t are the spatial and temporal coor- 
dinates in the comoving frame. $ is the electrostatic potential 
in the plasma, which is self-consistently obtained by solving 
Pois son's equation lTT2l 
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where 7 2 = 1 + \a\ 2 /2 as before. All quantities in eqs. ^ 
and ^ are normalised to the characteristic plasma parameters, 
i.e. space to c/uj p and time to 1/oj p . The intensity distribution 
a 2 and the refractive index rj obtained from the numerical in- 
tegration of eqs. ^ and ^ is shown in Fig. [4] 

The model shows that, at low densities, pulse compression 
starts from the back of the pulse, as observed experimentally. 
At I = 8 mm, the pulse has compressed to r w 21 fs, in 
good agreement with our measured compression. Laser en- 
ergy which sits in following buckets is trapped to form trailing 
pulses. By contrast the front of the pulse evolves only slowly. 
Since for relativistic pulses it is the front edge of the pulse 
which determines the wakefield amplitude, this also evolves 
slowly. Hence in this regime (low n e , high a), pulse compres- 
sion can extend the growth of the plasma wave to distances 
well beyond a naive depletion length obtained by considering 
only energy transfer to the wake. Furthermore here, depletion 
certainly cannot be modelled by pulse front erosion ifTTlfTTl . 

At these densities, there is little evidence for the explosive 
increase in a 2 predicted for propagation longer than the non- 
linear modification time 0, which can be stated simply as, 
I > cr n i « (cro/\a\ 2 )(n c /n e ) ifTTl . For n c /n e ~ 1000, and 
our initial pulse length, cr n i > 10 mm, and for I — 4 mm, 
this time is not reached for n e < 6 x 10 18 cm -3 , provided 
in both cases that \a\ 2 is not too much greater than 1. This 
supports the measurement shown in figure 3, that for our con- 
ditions, only a fraction of the laser energy is captured in the 
wakefield driving filament beyond zr. This also explains the 
good agreement with the ID calculations and also explains 
why there is not sufficient plasma wave growth to observe 
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FIG. 4: (color online) Numerical solution of ID wave equation in 
quasistatic approximation, a 2 (blue filled), n e (black solid), and r\ 
(red) resulting from propagation of a pulse with ao ~ 1.5 at n e o = 
2.5 x 10 18 cm -3 ; (a) initially, (b) at / = 8 mm. The reduced ao was 
chosen to account for the initial energy loss in the guided mode as in 
Fig.0 



self-injected electron beams below this density. Pulse evolu- 
tion leading to a 2 amplification and the resultant wavebreak- 
ing of the plasma wave has been inferred from the proper- 
ties of generated electron beam at the high end of this density 
range but over longer interaction lengths |2|. Though chal- 
lenging, future studies into this non-linear regime of pulse 
compression, which may produce extremely short laser pulses 
of higher power, may prove to be particularly rewarding. 
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